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We create a gray optical lattice structure using a blue de- 
tuned laser field coupling an atomic ground state of angular 
momentum J simultaneously to two excited states with angu- 
lar momenta J and J — 1. The atoms are cooled and trapped 
at locations of purely circular polarization. The cooling pro- 
cess efficiently accumulates almost half of the atomic popula- 
tion in the lowest energy band which is only weakly coupled 
to the light field. Very low kinetic temperatures are obtained 
by adiabatically reducing the optical potential. The dynamics 
of this process is analysed using a full quantum Monte Carlo 
simulation. The calculations explicitly show the mapping of 
the band populations on the corresponding momentum inter- 
vals of the free atom. In an experiment with subrecoil mo- 
mentum resolution we measure the band populations and find 
excellent absolut agreement with the theoretical calculations. 

PACS: 32.80.Pj, 33.80. Ps, 42.50.Vk 



I. INTRODUCTION 

Remarkable progress has been made in cooling and 
trapping of neutral atoms by laser fields. An outstand- 
ing success was the demonstration of Bose-Einstein con- 
densation |l| 0], which was achieved by combining laser 
cooling techniques with magnetic trapping and evapora- 
tive cooling. The quest to observe quantum many body 
effects by manipulating atoms merely using laser light 
has spurred researchers to invent new cooling techniques. 
Extremely narrow momentum distributions have been 
demonstrated by velocity-selective coherent population 
trapping Raman cooling [Q and adiabatic cooling 

|§|-|lol|. Subrecoil cooling of atoms in mesoscopic traps 
has been proposed (e.g. in Ref. p] , ^2[ ) and recently suc- 
cessfully demonstrated yielding a dramatic increase in 
the phase space density of the atomic vapor [ ^3j . 

Atoms interacting with spatially periodic light-induced 
potentials can be accumulated in an array of microscopic 
traps. The quantum motion of atoms in these optical 
lattices has been subject of extensive research Q]. The 
fascinating perspective that quantum statistical effects 
might become observable in optical lattices has stimu- 
lated the search for new types of optical lattices. Very 
promising schemes aiming at high atomic densities are 
gray optical lattices in which the trapped atoms are al- 



most decoupled from the light field [f5|. This reduces 
the density limitation by light induced atom-atom inter- 
actions. A one dimensional gray lattice using a laser field 
in combination with a magnetic field has recently been 
proposed |ll| and theoretically investigated [fill. Exper- 
imental demonstrations of magnetic field induced dark 
lattices in two and three dimensions are described in Ref. 

0- 

In this work we present an alternative and very simple 
way to efficiently accumulate atoms in a gray optical lat- 
tice structure. Our scheme merely uses a single frequency 
laser field and can be adapted to many alkali atoms. 
The adiabatic release of atoms trapped in the optical lat- 
tice is studied both theoretically and experimentally. We 
use a full quantum Monte Carlo wave function simulation 
to follow the atomic evolution during the lowering of the 
potential. It is shown that the band populations are in- 
deed mapped on the corresponding momentum intervals 
of the free atom, as was suggested by Kastberg et al. || 
for adiabatic cooling in a bright optical lattice. This pop- 
ulation mapping is accurate if optical pumping between 
different bands is negligible during the adiabatic release. 
In an experiment with subrecoil momentum resolution 
we demonstrate the population measurement of individ- 
ual bands using our gray optical lattice. For a wide range 
of parameters we find excellent agreement between the- 
oretically calculated band populations and the measured 
values. This establishes population mapping by adia- 
batic release as an experimental tool to investigate the 
interaction of atoms with a spatially periodic potential. 

We give a theoretical description of the gray optical 
lattice in section II and calculate its eigenfunctions in 
section III. Numerical simulations of the stationary en- 
ergy and position distribution using rate equations and 
quantum Monte Carlo wave function techniques are pre- 
sented in section IV. The adiabatic release is theoretically 
investigated in section V and experimentally in section 
VI. Section VII gives an outlook to experiments with high 
atomic densities. 



II. OPTICAL LATTICE CONFIGURATION 

We consider a one dimensional laser field acting simul- 
taneously on two atomic transitions coupling an atomic 
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ground state manifold with total angular momentum 
J g = J to two different excited levels with angular mo- 
menta J ei = J and J e2 = J — 1. The field consists of 
two counterpropagating waves with mutally orthogonaly 
linear polarization (lin_Llin). Fig. |l| illustrates the spa- 
tially varying polarization of the standig wave. The light 
field shall be detuned to the blue of both transitions and 
the resulting optical couplings shall have the same order 
of magnitude |y|. This can be realized on the D\ line 
of the alkali atoms Rb, Na and Li. The two excited lev- 
els are formed by hyperfine manifolds with the splitting 
Ahf '■= (Ej ei — Ej e2 )/h. An essential part of the system 
is the off-resonant coupling of the J e2 manifold with the 
detuning Aj + Ahf- 

The interaction of the oscillating atomic dipolc with 
the standing wave causes spatially periodic light shifts in 
the atomic ground state manifold. In regions of purely 
<j + (cr~) polarized light the atoms are optically pumped 
into the m= J g (m = — J g ) ground state, which is decou- 
pled from the light field and experiences no light shift 
or optical excitation. At locations of linearely polarized 
light all ground state sublevels are coupled to the ex- 
cited states and are shifted towards higher energies. In 
this semiclassical picture we expect that the atoms are 
cooled by a Sisyphus mechanism and accumulated in dark 
states at locations of pure cx-polarization. In a picture 
which treats the atomic motion quantum mechanically 
the atomic wavefunction always has a finite spatial ex- 
tend and can not be completely decoupled from the light 
field. We therefore expect the formation of a gray opti- 
cal lattice. Our situation is qualitatively similar to that 
of magnetic field induced dark optical lattices, where lo- 
calized dark states are created by combining a standing 
wave with a magnetic field [ pi|Jr7| ]. In the following we 
outline our calculations using the example of a J = 1 
ground state. This allows to demonstrate the essential 
physics using a minimal number of Zeeman sublevels. 

A model Hamiltonian for the J g = 1 — > J ei = 1, J e2 = 
system, which includes the kinetic energy is found in 
Ref. [p0[ . To simplify the following discussions, we as- 
sume low saturation, which allows to adiabatically elim- 
inate the excited atomic states Then the Hamilto- 
nian can be written as: 
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where M is the atomic mass, k the modulus of the light 
wave vector and \m) the m-th Zeeman substate of the 
ground state manifold. Uj ej — Aj ej sj ej /2 is the effective 
interaction potential (=light shift), where Aj ej — wi asor — 
uj Jej is the detuning and sj C] = Qj e /(2Aj Cj + 7^/2) is 
the saturation parameter for the transitions between the 
ground state manifold and the excited state manifolds. 
The optical pumping rates for the ground states are pro- 
portional to the parameter Tj ej — sj ej "/j ej , where jj e . 
are the lifetimes of the excited levels. 



III. OPTICAL POTENTIALS AND 
EIGENFUNCTIONS 

To calculate the adiabatic potentials we diagonalize 
the atom-field Hamiltonian Haf at each spatial point 
x separately. This yields the spatial dependence of the 
eigenvalues and eigenstates. In a semiclassical picture 
these eigenvalues amount to the light shifts experienced 
by an atom at rest in the corresponding states and acts 
as a potential for slowly moving atoms. Fig. |2| shows 
these adiabatic potentials for the case of a J g = 1 to 
J e . = 1,0 transition. Due to the threefold degeneracy 
of the unperturbed J g = 1 atomic ground state, we 
find three adiabatic potentials. The energetically low- 
est potential curve exhibits minima of zero light shift at 
locations of purely circular polarization. At the same 
spatial points x the maxima for the highest potential 
are found. The constant potential can be attributed 
to atoms in the m = state, which experience a con- 
stant light shift in space and hence feel no (semiclassi- 
cal) force. Atoms in state m = 1 are drawn towards re- 
gions of cr + -polarization, where they experience minimal 
(=zero) light shift. Atoms in the m = — 1 are repelled 
from this <r + -area (uppermost curve). The magnitude of 
both forces is comparable and an atom moving in the 
cr + -region spends most of its time in the weakly coupled 
m = 1 state. Therefore the semiclassical picture predicts 
a trapping force towards areas of circular polarization. In 
addition the corresponding optical pumping rates show 
the right spatial dependence to provide a Sisyphus cool- 
ing mechanism for moving atoms, as they are pumped 
into the locally less coupled (i.e. the energetically lower 
lying) states. 

We now include the motional degrees of freedom and 
calculate the eigenstates of the full Hamiltonian H. The 
diagonalization is performed numerically on a discrete 
spatial grid extending over several wells of the optical 
potential with periodic boundary conditions. The calcu- 
lated discrete eigenstates correspond to the energy bands 
of the full spatially periodic lattice structure. At energies 
above the well depth we find delocalized wavefunctions 
(unbound states) whereas at low energies the states are 
well localized (bound states). The atomic position dis- 
tribution \i/: g (x)\ 2 of the energetically lowest eigenstate 
is plotted with a dashed line in Fig. 0. The chosen ver- 
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tical offset is the ground state energy. As expected we 
find a strong localization of the wavefunction at the po- 
tential minimum of each well. Its width is a fraction of 
the optical wavelength and depends on the detuning and 
the intensity of the light field. The momentum spread 
associated with the finite size of the optical wells pre- 
vents the existence of eigenstates exactly decoupled from 
the light field. Nevertheless we find, that the localized 
states of lowest energy exhibit only a very small optical 
coupling to the two excited levels and hence they have 
correspondingly low light shifts and energies. 

The energy intervals between the lowest states are con- 
siderably larger than the recoil energy and reach up to 
a few hundred recoil shifts. Experimentally they deter- 
mine the Raman resonances found in fluorescence or in 
probe absorption spectra |l4] ]. With increasing energy 
the interval sizes decrease due to the anharmonicity of 
the potential. This leads to inhomogeneous broadening 
of the probe absorption resonances. 

The tunnel coupling between two equivalent states of 
neighbouring potential wells can be estimated from the 
eigenenergies of a two well calculation. The coupling 
leads to an energy splitting between the corresponding 
symmetrical and antisymmetrical states. For the param- 
eter regions and time scales we consider in the following, 
the tunnel coupling between the bound low energy states 
of the single wells is so small, that we can view the po- 
tential wells as independent. For situations with finite 
tunnel coupling it would be advantageous for calculations 
to use Bloch eigenf unctions of the optical lattice as the 
numerical basis 



IV. ENERGY AND POSITION DISTRIBUTIONS 

We now calculate the steady state distributions of 
atoms in the optical lattice. They allow to judge the effec- 
tiveness of the involved cooling and confinement mech- 
anisms on the basis of quantities as the mean energy, 
the position spread or the population of the energetically 
lowest bound states. 

The inclusion of the spontaneous scattering of photons 
in addition to the coherent atom-laser field dynamics im- 
plies a dynamical redistribution of the population among 
the various atomic states eventually leading towards a 
steady state |22| 20| , |TT[ ] . We calculate the steady state 
population distribution in a similiar approach as demon- 
strated in Refs. jT^|l^,|2(J using two methods: a rate 
equation apporach based on the Raman transition matrix 
elments between the eigenstates (2^Jl^] and a Quantum 
Monte Carlo wave function simulation technique (QM- 
CWFS) relying on a Bloch-state expansion of the atomic 
wavefunctions [^0|,^l). This allows to realistically treat 
periodic spatial geometries still maintaining a numeri- 
cally tractable grid size. We expect the rate equations to 
give fairly accurate results for the bound states, where 
the tunnel coupling is small, and less accuracy for states 



with energies near and above the potential well depth. 

The stationary position and momentum distributions 
obtained by the QMCWFS are shown in Fig. § The 
upper plot shows the total position distribution (solid 
curve) as well as the contributions of the various Zceman 
sublevels (dashed curves). As expected the m = state 
is almost not populated and shows no spatial variation. 
The localization towards points of purely circular polar- 
ization is state selective. We find better localization with 
increasing field strength. Due to the strong spatial con- 
finement to roughly Ax = A/8, we find a correspondingly 
large width of the momentum distribution of 5p > 8Hk, 
which is shown in the lower plot of Fig. ||. The obtained 
results are in excellent agreement with the results of the 
rate equation model as shown below. 

The optical pumping rates scale with the parameters 
Tj ej , which depend on the lifetimes of the excited states 
and on the detuning and the intensity of the laser field. 
Their relative magnitude can be tailored to a large extend 
by a suitable choice of these two laser parameters. The 
possible ranges of obtainable light shifts cx Uj e and opti- 
cal pumping rates cx T j e . may be limited by the available 
laser intensity and by the magnitude of the atomic hy- 
perfine splitting. An unfavourable choice of parameters 
can lead to additional unwanted off-resonant couplings 
to other atomic hyperfine levels. 

The influence of the lightshift Uj ej and optical pump- 
ing rate parameters Tj ej is demonstrated in Fig. ^| for 
the J g = 2 — > J e . —2,1 transition of the 87 Rb Di-line. 
We chose this transition as an example, because the hy- 
perfine splitting allows a wide range of this parameters. 
The figure shows the occupation probabilities for three 
sets of parameters. The data points marked with '*' cor- 
respond to a situation, where the coupling to the excited 
state manifold J ei = 2 is near resonant and the cou- 
pling to the J e2 = 1 manifold is off resonant. The data 
points marked with '+' correspond to a situation with 
couplings of equal strengths. The highest ground state 
occupation probability of 43% was found in the first case 
with Ui = 500 E R , U 2 = 1000 E R and ri/r 2 = 0.02. 
In all cases a large fraction of the atomic population is 
concentrated in the lowest few (gray) states. This cor- 
responds to strong spatial localization of the atoms, a 
result, which we obtain for the J g = 1 — > J ej = 1,0 
transition also by the QMCWFS. The computation time 
required for QMCWFS for transitions with J g > 1 is so 
high, that we selectively performed Monte Carlo simula- 
tions to verify the rate equation results. For transitions 
with higher angular momentum J g the population of the 
ground state (lowest band) is higher, because the rel- 
ative magnitude of the Clebsch-Gordan coefficients for 
the states with m — ±J g and m = ±(J g — 1) changes. 

The energy distributions in the example strongly de- 
viate from thermal distributions of the same mean en- 
ergy. Therefore we do not have a thermal equilibrium 
state to which one could consistently attribute a tem- 
perature. The ground state population for the case with 
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U 2 = 1000 E R and r x /r 2 = 0.02 is w 20% higher than 
for a thermal distribution of same mean energy (E) « 
130 En. So the temperature deduced from the relative 
populations of the lowest two levels kel^ = 76 Er differs 
by w 25% from the corresponding value keT^ = 96 Er 
obtained from levels n = 3 and n = 4. The disagreement 
becomes more significant for higher n. 

The mean energy (E) of the atoms, which is e.g. of 
central importance for loading of a purely magnetic trap, 
and the stationary value of the atomic population Pi of 
the lowest energy band are key quantities with respect 
to a possible direct observation of quantum statistical 
effects. The dependence of (E), Pi and P2 on the light 
shift U2 for fixed ratios ri/I^ = 0.1 and Ui = U 2 j2> is 
shown in Fig. ||. This case corresponds to a laser tuned 
far to the blue of the J g — 2 — > J e . =2,1 transitions. A 
probability of Pi = 45% to find the atom in the lowest 
energy eigenstate is achieved for parameters well in the 
reach of experimental capabilities. 

V. ADIABATIC RELEASE 

It has been suggested in recent experimental work 
P,[l0| that the populations of an optical lattice can be 
directly experimentally measured by an adiabatic release 
of the atoms from the lattice and a subsequent measure- 
ment of the resulting atomic velocity distribution. If the 
release is fully adibatic, the atoms from the lowest (first) 
band will be mapped exactly to the momentum inter- 
val between —hk and +hk and the second band will be 
mapped to the intervals —2hk to — hk and Tik to 2hk. The 
n-th band will be mapped on the —nhk to — (n — l)hk 
and the (n — l)hk to nhk intervals. Nonadiabaticity and 
incoherent redistribution of the atoms during the release 
will alter the mapping and the assignment of the band 
populations of the lattice and the momentum intervals of 
the free atoms will be less accurate. We have performed 
a QMCWFS with a time dependent laser intensity to 
quantitatively verify the mapping between the stationary 
population distribution of the lattice and the momentum 
distribution of the free atom. All the effects of nonadia- 
baticity and incoherent spontaneous redistribution of the 
atoms during the release are fully accounted for in this 
model. 

The time evolution of the atomic momentum distri- 
bution is shown in Fig. |^ for an gradual turnoff of the 
lattice field. For t — we start with the steady state 
distribution, which we calculated by the rate equation 
approach disussed above. The atoms are localized in the 
lattice and correspondingly the momentum distribution 
is broad. When the lattice field intensity is reduced, the 
momentum distribution becomes narrower due to adia- 
batic cooling. For the given example the optical potential 
varies as Uj ej (t) = Uj ej (0)exp(—(t/T) 2 ) as it would oc- 
cur if an atom leaves a Gaussian beam transversly at a 
constant velocity. The chosen time scale of the turnoff 



t = V2tr (where tr — 2nh/Eji is the recoil time) is 
long compared to the oscillation period of atoms in the 
lattice wells but short compared to the lifetime of the low- 
est dark levels. The lightshift parameters Ui = 200 -Er, 
Uq = 100 En and the optical pumping rates Ti — 3 Er/U 
and To = 1 Er/Tl are the same as in Fig. ||. 

The narrowing of the momentum distribution contin- 
ues for f>r down to a width of roughly one recoil mo- 
mentum. From the resulting momentum distribution we 
calculate the probability Wn to find an atom in the in- 
terval between -Nhk and +Nhk. We find W x = 37, 9%, 
W 2 = 56,4% and W 3 = 68,1%. This agrees to within 
less than 1% with the corresponding initial (steady state) 
populations of the lattice Pi = 37, 4%, Pi + P 2 = 56, 3% 
and Pi + P% + P3 = 67, 9%. A comparison for higher 
momenta is limited by the accurracy of the rate equation 
approach. 

VI. EXPERIMENT 

To demonstrate the experimental feasibility of the pro- 
posed method to measure the band populations of the 
optical lattice, we performed an experiment with the ap- 
paratus described in Refs. ]^Jl0[|. Our gray lattice scheme 
is realized on the J g = 3 — ► J ej =3,2 transitions of the 
85 Rb Z?i-line (where J is the total angular momentum of 
the atom). 

A pulsed beam of cold Rubidium atoms is directed 
vertically downwards and crosses a standing wave field 
(lattice field), which induces the optical potentials of the 
dark lattice. The atoms are cooled into the lattice sites 
and are then gradually released from the optical potential 
due to the Gaussian shape of the lattice field. Below the 
lattice field the momentum distribution of the atoms is 
measured with a resolution of one third of the photon 
recoil. 

The lattice field is induced by a standing wave oriented 
along the x-axis and has a frequency tuned A3 = 2673 
to the blue of the J g — 3 — > J ei =3 transition. The 
hyperfine-splitting between the two excited states of the 
-Di-line is Arr = 65 73 and the detuning of the sec- 
ond excited states is A2 = 91 73. The incoming beam 
of the standing wave is linearly polarized along the z- 
axis and the reflected beam is polarized along the y-axis. 
The Gaussian waists of the beams are w z — 1.35 mm in 
z-direction and w y = 0.48 mm in y-direction. This corre- 
sponds to a 0.4 ms time of flight of the atoms (3.2 m/s) 
through the waist w z . The region of the lattice field 
is shielded against stray magnetic fields to well below 
0.5 mG. A second standing wave overlaps the lattice field. 
It is tuned to the J = 2 — * J = 3 transition of the 85 Rb 
D2-line and optically pumps the atoms into the J = 3 
groundstate. 

To determine the atomic momentum distribution a 
pinhole, 75 urn in diameter, is placed 5 mm below the 
standing wave axis and the spatial distribution of those 
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atoms passing through the pinhole is measured 9.6 cm 
further down by imaging the fluorescence in a sheet of 
light. The sheet of light is formed by a standing wave 
which is resonant with the J = 3 — > J = 4 closed cy- 
cle transition of the D2-line. For each set of parameters 
we accumulate 200 single shot images and subtract the 
separately measured background. To obtain a one dimen- 
sional momentum distribution in x-direction we integrate 
the two dimensional distributions along the y-axis. 

Fig. Q shows an experimentally measured momentum 
distribution (dotted curve) in comparison with the cor- 
responding distribution obtained from the QMCWFS 
(solid curve). The initial distribution for the Monte Carlo 
simulation was the steady state distribution of the optical 
lattice in the center of the Gaussian beam. The calcula- 
tion was performed as explained in section ^ but using 
a finer momentum grid. The plotted curves are in good 
agreement except for a small fraction of atoms. In the ex- 
periment more atoms are found at higher velocities and 
less at low velocities. This can be related to the finite 
interaction time. The atoms in the experiment do not 
reach the steady state distribution in the center of the 
Gaussian beam and some fast atoms are not yet cooled 
into the lattice wells. 

For the experimental data we count the number of 
atoms detected in the momentum intervals —hk to +hk 
and —2hk to +2hk corresponding to the populations in 
the lowest energy band and the two lowest energy bands, 
respectively. These experimentally obtained populations 
are plotted in Fig. || with data points versus the lightshift 
U3 on the J g = 3 — > J ei =3 transition in the center of 
the Gaussian beam. The data points were recorded for 
several intensities and fixed detuning. The solid line rep- 
resents the (steady state) band populations in the lattice 
calculated for the center-of-beam parameters using the 
rate equation approach. The experimentally measured 
populations and the calculated steady state populations 
agree within 5% over the full investigated range of pa- 
rameters. This is remarkable, because the calculation 
was based only on the detunings and the intensities in the 
center of the Gaussian beam and the comparison involves 
no fit parameter. The small deviations for high lightshift 
parameters U3 towards higher ground state population 
for the experimental values can be attributed to a small 
but finite spontaneous emission probability during the 
release of the atoms, which especially effects the energet- 
ically higher lying less dark states and which transfers 
additional population to the ground state. 

VII. OUTLOOK 



well of an optical potential and that laser-like sources for 
matter waves might become feasible. In our theoretical 
work we have not included any atom-atom interaction. 
Thus we can only speculate that the low photon scatter- 
ing rate for atoms trapped in a gray optical lattice leads 
to much higher achievable atomic densities than for the 
case of bright optical lattices. This optimism is based on 
the assumption that atom-atom perturbations by dipole- 
dipole interaction and reabsorption of scattered photons 
are much reduced for atoms trapped in a dark state. The 
atomic densities necessary for the observation of quantum 
statistical effects are considerably lower in steep poten- 
tials compared to the cases of wide traps or free particles. 
The optical potentials of gray optical lattices can have an 
energy spacing much larger than the recoil energy and 
the trapped atoms can be cooled to mean energies of the 
same order of magnitude. We reached in one dimension 
a ground state occupation probability of ss50%, so that 
- extended to three dimensions — already two atoms 
located in the same potential well are sufficient for quan- 
tum statistical effects to become relevant. This might be 
observable even at average filling factors below one atom 
per lattice site. 

In conclusion, we have theoretically and experimen- 
tally studied a gray optical lattice structure which com- 
bines a low photon scattering rate with a high population 
in the lowest energy band. The lattice is formed by cou- 
pling an atomic ground state to two excited states. The 
atoms are trapped at locations of purely circular polar- 
isation which allows an extension of the scheme to two 
and three dimensions using the same field configurations 
as for bright optical lattices. We have numerically simu- 
lated the dynamics of atoms adiabatically released from 
the optical potential and the mapping of the band pop- 
ulations on the corresponding momentum intervals. The 
quantitative agreement with the band populations mea- 
sured in the experiment shows that adiabatic release is 
a promising tool to study the density dependence of the 
band populations in an optical lattice. 
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An interesting field for future experimental and theo- 
retical work is the interaction of a high density atomic 
sample |23| with a periodic optical potential p4| . It 
|2q^q| that quantum statistical ef- 



has been suggested 
fects may tend to cluster bosonic atoms within a single 
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FIG. 1. Schematic setup of the gray optical lattice. 
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FIG. 2. Adiabatic potentials Vdark, VbrigM and Vm=o (solid 
curves) and the atomic position distribution of the lowest en- 
ergy eigenfunction ip g shifted by the corresponding eigenen- 
ergy (dashed line) for the J g = 1 — + J ej = 1,0 system for 
Ui = 400 E R and U = 400 E R . 
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FIG. 3. Position and momentum 

distribution for Ui = 200 E R , Uo = 100 E R , I\ = 3E R /h 
and To = lE R /h for the J g = 1 — » J e . =1,0 system. The 
uppermost lines give the total distributions, while the lower 
lines give the partial contributions of the various Zeeman sub- 
states. The dotted curves correspond to the m = substate. 
The momentum distribution is normalized and given in units 
of (hk) -1 , the position distribution is periodic and given in 
arbitrary units. 



]2Q Energy and Population of the lowest two bands 
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FIG. 5. Ground P\ and first excited band P2 pop- 
ulations (short dashed curves), the corresponding energies 
Ei, E2 and the mean energy {E) (dot dashed line) for the 
J g = 2 — > J e . —2,1 system as a function of the laser intensity 
oc U2 for fixed ratios of the lightshift potentials U1/U2 = 1/3 
and of the optical pumping rates Fi/r2 = 0.1. 



Stationary population distribution 

0.45 



0,40 






1 000 E„ 


r,/r 2 =o.02 


0,35 




U,=500E„ +U 2 = 


1 000 E B 


r,/r 2 =i,o 









100 E R 


ryr,=o.i 


0,30 










0,25 










0,20 


* 










0,15 










0,10 











0,05 
0,00 




t 

* 

* * 

i i 8 


S 


* * © 




3 


5 

Band index 


10 


1 



FIG. 4. Population of the low energy eigenstates for the 
J g = 2 — > J ej = 2, 1 system for various lightshifts U2 and 
ratios of the optical pumping rates T1/T2 



Atomic momentum distribution during expansion 
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FIG. 6. Normalized atomic momentum distribution dur- 
ing various stages of the field turnoff t = 0, t = 2 r R and 
t — 4 T R for a J g = 1 — > J ej =1,0 system and for the same 
parameters as in Fig. [| The optical potential decreases as 
Uj ej {t) = Uj ej (0)exp( — ^(t/TR) 2 ). The units on the vertical 
axis are (Tik)" 1 . 
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FIG. 7. Normalized momentum distribution of the atoms 
after their adiabatic release from the lattice. The QMCWF 
simulation corresponds to the solid line, while the dashed line 
gives the corresponding experimental result for A3 = 26 73 
and U3 — 73/4. The units on the vertical axis are (Tifc) -1 . 
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FIG. 8. Population of the lowest (lower curve) and lowest 
two (upper curve) bands of the optical lattice as a function 
of the lightshift as obtained by quantum rate equations (solid 
lines) in comparison to the experimentally measured results 
(data points) for the J g = 3 — > J ej — 3, 2 transition of the 
8j Rb Di-line. The detuning was chosen to be A3 = 2673 to 
the blue of the J g = 3 — > J ei =3 transition. 
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